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We performed the calculation of the full O(a) corrections to e + e~ — > [lTvud with the help of the GRACE/1-L0QP 
system. We discuss how a finite decay width introduces a serious gauge invariance breaking, particularly for 
infrared 5-point functions. This is related to the way the reduction of those functions is performed and to the 
treatment of the width in the reduction. 



1. Introduction 

The first calculation of the full 0(a) correc- 
tions to e + e~ — * uVH was presented at the RAD- 
COR/L&L 2002 workshop Q. Since then, sev- 
eral authors published the radiative corrections 
to several important 2 — > 3 processes: e + e~ — * 
i/Fff |2l8l4lfi| , e+e~ -> UH \(il7im\ . e+e" -> 
ZHH [HSn\ . 77 -» ttFUa, e+e" -» e+e~H^ 
and e + e~ — > ^F^|13|. Now, full EW 1-loop calcu- 
lations are well under control for 2^3 processes 
in the SM. In this paper, we discuss the case for 
a 2 — > 4 process. 

As a first trial we take a typical LEP-2 process 
e + e _ — > [irVud. Although a status report on a 
similar attempt has appeared^!], there still is no 
complete result. 

2. Motivation 

For LEP2 experiments, the Double Pole 
Approximation(DPA) Ijj and the fermion loop 
scheme were used to predict the cross sections 
for e + e _ — > 4-fermions. They have the following 
features:(l) Gauge invariance is guaranteed. (2) 
It was sensible to split all 10 tree diagrams(CClO) 
into the doubly resonant ones(CC03) and others. 
The constant width is introduced in a naive way. 
In the energy region at the future linear-collider, 



however, non-CC03 diagrams are not negligible. 
For example, at s/s = 500 GeV, the tree level 
cross sections for CC03 and CC10 are 213.56 fb 
and 222.39 fb, respectively. The difference be- 
tween them reaches 4%. Therefore, the size of 
the radiative corrections to CC10 should be care- 
fully estimated at the TeV energy region. 

3. Structure of the calculations 

GRACE/l-LOOPUm has been used for the calcu- 
lation which proceeds through the following steps: 
1) Evaluation of the numerators, for which the 
symbolic manipulation system is used. In or- 
der to shorten the size of the matrix elements, 
the fermion masses m e , m M , m u , md are neglected 
when integrating over phase space. 2) Evaluation 
of the loop integrals. After 5- and 6-point inte- 
grals are reduced to 4-point functions, FF|17j and 
other analytic formulas are invoked for 4-, 3-and 
2-point integrals. All masses are kept in the loop 
integrals. 3) Construction of kinematics. Here 
masses are also kept exactly. 

4. The reduction algorithm of iV(> 5) point 
functions 

The amplitude with 5-point function 18 3 10 is 

T (5) f d ' 11 hjjj "' lp m 

J (2n)"i D D! ■ ■ ■ Di ' { ' 



1 



2 



7 

-wwvwvy 






wwwy — ' 


e 





Figure 1. Example of unstable 4-point functions 
obtained in the reduction of a 6-point function 



where D — l 2 — m 2 , and D t = (I + s^) 2 — m h * = 
1, • ■ • , 4 and Sj is a combination of external mo- 
menta. 

Multiplying the identities 

4 

i,i'=l 

with the loop momentum I, 
i" = ^ «f (A - (i • Sj), 



(2) 



where Ay — Si ■ Sj and Ai = sf — mf, we can 
reduce the numerator to 

ly.lv ■■■lp = ^^2(A- l ) ij (D i -D -Ai)s l ;i v ---l p . (3) 



This method is also applicable for the reduction 
of a 6-point function to a sum of 4-point func- 
tions. There is also a standard reduction method 
where the following identity is used; 



i 2 =J2v- *<w •«*■)• 



(4) 
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Because this method not only raises the rank 
of the momentum tensor in the numerator but 
makes the resultant source codes very lengthy 
compared with EqEJ GRACE/1-L00P uses this 
identity only as an option. As mentioned the sys- 
tem employs the FF package for the evaluation of 



4-point functions. It happens, however, that this 
package leads to numerical instabilities or incon- 
sistency in some cases having to do with some 
infra-red boxes. This occurs for instance when 
an internal massless particle is involved in some 
non-IR boxes like those obtained from a 6-point 
function as shown in Fig.l. Some infra-red boxes 
also need to be regulated by the introduction of a 
width for the particle(W or Z) circulating in the 
loop. We implement a constant fixed width for 
such cases. For all these particular cases we use 
special in-house routines. For the computations 
we performed till now these routines have worked 
quite well for any 2^2 and all 2 — > 3 one-loop 
processes. However, for the problem at hand, we 
have met a difficulty in carrying out the evalua- 
tion of a 5-point function when the W width is 
required. 



5. 1-loop diagrams for e + e — > /j, vud 

We introduce two sets of diagrams, the full 
set and the production set. In the standard 
model and within the non-linear gauge fixing 
conditions j one counts 44 tree- level and 6094 
1-loop diagrams. This defines the full set. 

When the masses of the electron, /i, u- and d- 
quark are ignored in the numerator of the ma- 
trix elements, we obtain the production set which 
consists of 10 tree diagrams(CClO) and 668 1- 
loop diagrams (with 88 pentagon and 60 hexagon 
graphs). 

6. Check of codes 

The following set of input parameters are used 
in the calculation: Throughout this paper, the 
results are expressed in terms of the fine struc- 
ture constant in the Thomson limit = 
137.0359895 and the Z mass M z = 91.1876GeV. 
The on-shell renormalization scheme uses M\y as 
an input parameter. Nevertheless the numerical 
value of My/ is derived through Ar(201- Mw thus 
changes as a function of My. We take m u = 
m d = 63 MeV, m s = 92 MeV, m c = 1.5 GeV, 
TO/,=4.7GeV. This set of masses gives a "pertur- 
bative" value of a(Mz) compatible with the cur- 
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Table 1 



Stability of the matrix element squared, Eq. 5, against the changes of Cuv and A. 



graph 


Cuv 


A = 10- 18 GeV 


A = 10-" 1 GeV 


full 
full 

prod. 

prod. 




100 


100 


-0.85388129300841993023002220162E-03 
-0.85388129300841993023002220166E-03 
-0.85388086448063959950419248157E-03 
-0.85388086452253511266116024563E-03 


-0.85388129300794622839025322460E-03 
-0.85388129300794622839025322465E-03 
-0.85388086990548652947683563335E-03 
-0.85388086994738204263380339741E-03 



rent value derived experimentally. We also take 
mt op = 180GeV. With these values, we find Mw 
=80.4163 and r z =2.4952GeV for Af ff =12QGeV. 
IV=2.118GeV is taken from PDG. 

The results of the calculation are checked by 
performing four kinds of tests at few points in 
phase space. For this we worked with the full 
set of diagrams in quadruple precision. The first 
check is the ultraviolet finiteness test. The regula- 
tor constant Cuv = l/ e — 1e + log An, n = 4 — 2e 
is kept in the matrix elements. In varying this 
parameter Cuv, we found that the result is sta- 
ble with an accuracy reaching 29 digits. In- 
frared(IR) finiteness is also checked by introduc- 
ing a fictitious photon mass A, treated as an in- 
put parameter in the code. The sum of loop and 
brcmsstrahlung contributions stays stable against 
variations in A with an accuracy of 12 digits. Ta- 
ble 1 shows the stability for both cases of full 
set and production se£(prod.). For the latter the 
accuracy worsens to order 0(m?/s) because the 
masses of the fermions are neglected in the nu- 
merator. 

The third check relates to the independence 
on the parameter k c which is a soft photon 
cut parameter that separates soft photon ra- 
diation(analytical formula) and the hard pho- 
ton performed by the Monte-Carlo integration. 
Gauge parameter independence of the result is 
performed as the last check through a set of five 
gauge fixing parameters. For the latter a gen- 
eralized non-linear gauge fixing conditional nas 
been chosen. In total the amplitude contains 5 
free parameters: 5,/3, S,e and k. In this check 
we set Tw = Tz = and keep all masses in the 
numerator. The whole matrix element shows no 
dependence on any of the gauge parameters in 
more than 21 digits. 
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Figure 2. An IR 5-point function and a tree dia- 
gram. 

7. Status of calculation 

All the calculations were done in quadruple pre- 
cision. When we sampled points in phase space 
to carry out the MC integration, we encountered 
an unstability with IR 5-point diagrams. In Fig. 2 
we show such an example with a product of a 
1-loop IR diagram and a tree diagram. It is to 
be noted that the 1-loop graph is obtained from 
the tree by adding a photon that joins e~ and 
u. A 5-point diagram is written as a sum of five 
4-point integrals. For the product of the two dia- 
grams (the labelling is assigned automatically by 
the system), we then write 



[IR-bop]^ = 2^245^ 

T (A)IR j(A)IR r (4) j{A) r (4) 



(5) 



which contains two infrared divergent boxes. Cur- 
rently the widths Tyv and T z are retained only 
in the IR part of the box integrals in such a way 
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Table 2 

Cancellation between each J^-point function in Eq.(8) at a phase point with the matrix elements of Fig. 2 



box 


IV = o 

(real, imag) 


IV 7^0 
(real, imag) 


1 

2 
3 
4 
5 


( -17.69217505, 66.88369265) 
(-929.96695924, 1372.81890328) 
( 782.70500402, -5690.37576950) 
( 410.83071602, 4557.14502347) 
(-245.31384970, -301.36098220) 


( -4.22744975, 2.9465712570) 
( -120.36833308, 39.2204483424) 
( 321.80994609, -283.271079293) 
( -199.63435655, 275.873356743) 
( 0.65534291, -27.683645798) 


sum 


( 0.56273604, 5.11086770) 


( -1.76485038, 7.085651250) 



and also in the reduction formulas for 4-points. 
Without the width, 3 digits cancellation is ob- 
served among the boxes, as shown in Table 2 at 
some phase space point near the W^-pole, 

(q 2 - + iT w M w ) ud 

= -169.3611 + 167.1833i (GeV 2 ) (7) 

On the other hand, when IV 7^ 0, Table 2 shows 
less cancellation. This led to a bad convergence 
of the phase space integration as we will clearly 
demonstrate below. Though the mechanism for 
this cancellation with zero width is not fully un- 
derstood now, it is highly probable that an iden- 
tity including 

{(l.p 1 )e(p 2 ,P3,P4,P5) ~ (l-P2)e(pi,p 3 ,P4,Pb) 
+ (l-P3)e(Pl,P2,P4,P5) - (l-P4)e(pi,P2,P3,P5) 

+ {l.p 5 )e(p 1 ,p 2 ,P3,P4)} x (l.pe) = 

causes the trouble. This identity relates the five 
box integrals among each other. No such phe- 
nomenon took place for 6-point diagrams nor non- 
Ill 5-points. Also we did not encounter such dif- 
ficulty in all 2 — * 3 processes we studied till now. 
The number of independent fermion lines may be 
related to this undesirable situation. As men- 
tioned above, GRACE/ 1 -LOOP adopts Eq|3|for the 
reduction. We found, however, that when we used 
EqUjno such problem happened. 

8. Test run of integration 

In order to confirm that an improper treat- 
ment of the width in the infrared divergent 5- 
point function is the cause of bad convergence, 
we temporarily took the following ad hoc regular- 



ization in the loop amplitude 

- (5) - 25R ( T IRW> T tree) 

, 9S p - M W ( rp \ T t 

[qi - l J IR-W 5 Viv=o tree 

(8) 

and integrated all the diagrams in the phase 
space. The integration is carried out by BASES. 
Those diagrams which contribute less than 0.01 fb 
were omitted to leave 361 1-loop diagrams. Hard 
photon emission is also included. 

At = 500 GeV, a tree = 222.39 ±0.01fb. Ta- 
ble 3 shows the results of the integrated cross sec- 
tion for each set of n-point diagrams with 50,000 
sampling points, together with their MC error. 
The column "original(IR)" shows the results of 
the diagrams related to the infrared divergence 
as extracted normally through Eqs. [SJ-6 or an 
equivalent form for the other n— point IR dia- 
grams. On the other hand "ad /ioc(IR)" is the 
result after making the modification in EqlHl for 
the same sets of diagrams as "original(IR)" . The 
last column contains those diagrams which do not 
have any infrared divergence. For the case of 5- 
point functions of the type "original(IR)" , the in- 
tegrated error is extremely huge due to the mis- 
match between its corresponding 4-point parts as 
mentioned earlier. After applying EqlHl there is 
a drastic improvement in the error and the cen- 
tral value of the integration is also shifted. Other 
contributions from n-point functions remain the 
same within the integration error if one applies 
the 11 ad hoc" prescription. Table 3 suggests that 
a proper treatment of the width in 5-point func- 
tions is crucial to improve the situation but one 
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still needs to find a proper justification for the 
factorization. 

Table 3 

Integrated cross sections in fb showing the contribu- 
tion of each set of n-point diagrams for the test run. 
The parentheses show the MC integration error. 
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9. Summary 

One-loop amplitudes of the full 6094 diagrams 
for e + e~ — > [irvud were generated by the 
GRACE/1-L00P system. Non-linear gauge invari- 
ance has shown the consistency of the full set of 
amplitudes and the system itself. A new reduc- 
tion algorithm from a 6-point function to 4-point 
function works well. A finite decay width brings 
a serious breaking of gauge invariance, particu- 
larly for 5-point infrared integrals. It is clear now 
that the radiative corrections to 2 — > 4 processes 
are calculable, though more improvements are in- 
evitable. 
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